%% term_U_V
function matrix_A = term_U_V(mesh, fespace)
A1 = term_P_Q(mesh, fespace);
matrix_A = blkdiag(A1,A1);
end